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We prove that the Mott insulating state is characterized by a divergence of the electron self 
energy at well-defined values of momenta in the first Brillouin zone. When particle-hole symmetry 
is present, the divergence obtains at the momenta of the Fermi surface for the corresponding non- 
interacting system. Such a divergence gives rise to a surface of zeros (the Luttinger surface) of 
the single-particle Green function and offers a single unifying principle of Mottness from which 
pseudogap phenomena, spectral weight transfer, and broad spectral features emerge in doped Mott 
insulators. We also show that only when particle-hole symmetry is present does the volume of 
the zero surface equal the particle density. We identify that the general breakdown of Luttinger's 
theorem in a Mott insulator arises from the breakdown of a perturbative expansion for the self 
energy in the single-particle Green function around the non- interacting limit. A modified version of 
Luttinger's theorem is derived for special cases. 



I. INTRODUCTION 

In the absence of disorder, electronic states in insula- 
tors fail to carry current either because the band is full (a 
band insulator) or a gap opens at the chemical potential. 
The latter is indicative of either a structural transition in 
which a partially filled band splits by doubling the unit 
cell or a Mott state in which strong on-site correlations 
create a charge gap in a half-filled band, as in transition 
metal oxides, the high-temperature copper-oxide super- 
conductors a case in point. While the origin of the Mott 
gap is clear, the mechanism for the bifurcation of the 
half-filled band into upper and lower Hubbard bands re- 
mains controversial^. We show here that the Mott gap 
leads to a divergence of the self-energy in momentum 
space, which in the case of particle-hole symmetry lies on 
the non-interacting Fermi surface. Such a divergence 
prevents the renormalized energy band from crossing the 
Fermi energy and hence is ultimately the mechanism by 
which Mott insulators insulate. Further, this divergence 
leads to a vanishing of the single-particle Green function 
thereby defining a surface of zeros in the first Brillouin 
zone. We demonstrate that the zeros account for nu- 
merous anomalous properties of the normal state of the 
cuprate superconductors. The zeros found here obtain 
fundamentally from strong coupling Mott physics'* and 
should be distinguished from those arising from weak- 
coupling symmetry-breaking instabilities^ of the Fermi 
surface. 

That the Fermi surface of the non-interacting system is 
in some way connected with strong-coupling Mott physics 
(that is, Mottness) has not been anticipated. For Fermi 
liquids, such a correspondence is natural. Fermi liquid 
theory^ rests on the seemingly simple conjecture that the 
number of bare electrons at a given chemical potential 
equals the number of Fermi excitations (quasiparticles) 
in the interacting system regardless of the strength of 



the interactions. Quasi-particle excitations are identified 
by simple poles of the single-particle, time-ordered Green 
function, G{p,e). Hence, for a Fermi liquid, the Landau 
conjecture, proven by Luttinger—, implies the existence 
of a surface (the Fermi surface) in momentum space at 
which 5RG'(p, e = cf) changes sign by passing through 
infinity. In systems lacking quasiparticles (no poles), 
such as insulators, the Landau correspondence between 
the particle density and quasiparticle excitations breaks 
down. In fact, as uncloaked recently^i^ii, Luttinger's the- 
orem is not necessarily invalidated when quasiparticles 
are absent. The suggestion^iii^ is that the particle den- 
sity 

N ^ f (Pp _ 

^ Jg{p,o)>o {^^r 

is properly defined through an integral in momentum 
space over a region where the single particle propagator 
is positive. An explicit assumption in Eq. ([T]) is that the 
imaginary part of the self-energy vanishes at the chem- 
ical potential^. Sign changes occur at poles or at zeros 
of the Green function. Should Eq. ^ hold, then the 
volume enclosed by the surfaces of zeros and poles yields 
the particle density. 

In this paper, we offer a criterion for the location of the 
surface of zeros and show explicitly that the volume en- 
closed by the zeros yields the electron density only when 
a Mott insulator has particle-hole symmetry. Our proof 
of the latter makes general the perturbative arguments 
made recently for a Kondo insulator—. In the absence of 
particle-hole symmetry, there is an additional contribu- 
tion to the electron density in Eq. ^ which arises from 
the breakdown of perturbation theory. We show explic- 
itly that the breakdown of Luttinger's theorem does not 
obtain from aT — regularization of divergent integrals 
as has been claimed^ii^ but rather arises anytime the 
self energy cannot be obtained perturbatively around the 
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non-interacting limit. Our results are in agreement with 
the mechanism proposed by Altshuler, et. al^i^ in the 
context of the breakdown of Luttinger's theorem in the 
presence of a spin-density wave. Finally, we demonstrate 
that models which project out the high-energy scale at 
half- filling lose the surface of zeros. As a consequence, 
the quasiparticle weight need not^^ be the same in the 
Hubbard and projected schemes such as the t-J model. 
The evolution of the surface of zeros in the pseudogap 
phase is also delineated. 

II. ZEROS IN A MOTT INSULATOR 

As our starting point, we consider the simplest model 
which captures the physics of Mott insulators, the Hub- 
bard model, 

iJ = - ^ Ujcl^Cja + UY^ n,^n,i - rij^, (2) 

in which electrons hopping on a lattice between neigh- 
bouring sites with amplitude = toy and chemical 
potential fi pay an energy cost U anytime they doubly 
occupy the same site. The operator Cicr (cjg.) annihilates 
(creates) an electron on site i with spin a and /x sets the 
chemical potential. The quantity of interest is the single- 
particle retarded Green function, G'^^{i,j,t,t') = —i9{t— 
t'){{cicr{t), Cja-it')}), in particular, its momentum and en- 
ergy Fourier transform, G™*(fc,e) = FTG'^\i, j,t,t'), 
where {a, 6} indicates the anticommutator of a and b and 
9{x) the Heaviside step function which is non-zero only if 
its argument is positive. The quantity that is directly ob- 
servable experimentally through ARPES is the spectral 
function, ^^(k, e) — — 3GJj°'(k, e)/7r. Summed over mo- 
mentum, the spectral function defines the single-particle 
density of states. The causal nature of the Green function 
permits it to be constructed entirely from its imaginary 
part 

Gr(k,c.)= Tdc^'AMI (3) 

through the standard Hilbert representation. For a Mott 
insulator, a gap of order U occurs in the spectral function. 
We will take the gap to have a width 2A centered about 
0. As we consider the general case in which symmetry- 
breaking plays no role in the gap, simply strong electron- 
correlations, our conclusions regarding the zeros are ap- 
plicable to the mechanism proposed by Mott. Within 
the gap, A(k, e) = 0. This is a necessary condition for 
any gap. Consequently, in the presence of a gap, the real 
part of the Green function evaluated at the Fermi energy 
reduces to 

ii;,(k, 0) = - /'^^ de'^^^^ - r (4) 

as pointed out by Dzyaloshinskii-. At half-filling, Mott 
insulators have half the spectral weight above the gap. 



Hence, it is, possible, in principle, that the real part of 
the Green function vanishes along some momentum sur- 
face. However, no criterion has been given for the loca- 
tion of such a surface. Indeed, the nature of the Mott 
transition has been studied extensively in d = oo us- 
ing dynamical mean-field theory (DMFT)i^. In DMFT 
{d — oo), the self-energy diverges at w = for all mo- 
menta. Hence, in d = oo, there is no Luttinger surface. 
For finite dimensional systems, it is not known what be- 
comes of the unphysical momentum-independent 1/lu di- 
vergence of the self energy in d = oo The new feature 
that this work brings into focus is the divergence of the 
self energy along a continuously connected momentum 
surface in the first Brillouin zone as the defining feature 
of a Mott insulator in finite dimensions. 



A. Particle-Hole Symmetry 

We now prove that when particle-hole symmetry is 
present, the spectral function is an even function of fre- 
quency at the non-interacting Fermi surface. As a conse- 
quence, Eq. (HI) is identically zero along that momentum 
surface. To proceed, we consider a general particle-hole 
transformation, 

Cia — > e'*^ '''cJo-> (5) 

of the electron annihilation operator. That the Hamil- 
tonian remain invariant under this transformation places 
constraints on both Q and the chemical potential. The 
Hubbard model with nearest-neighbour hopping remains 
invariant under Eq. ([5]) for Q = (tt, it) and /i = U/2 . 
The latter is the value of the chemical potential at half- 
filling, the Mott state. Transforming the operators in 
the Green function according to Eq. ([5]) and keeping the 
chemical potential fixed at /i = C//2 leads to the identity 

A^(k,uj) = A^{-k-Q + 2mT,-uj). (6) 

Hence, the spectral function is an even function of fre- 
quency for k = Q/2 -|- riTT. Consider one dimension and 
nearest-neighbour hopping. In this case, the symmetry 
points are ±7r/2, the Fermi points for the half-filled non- 
interacting band. In two dimensions, this proof is suffi- 
cient to establish the existence of only two points, not a 
surface of zeros. To determine the surface, we take advan- 
tage of an added symmetry in higher dimensions. For ex- 
ample, in two dimensions, we can interchange the canoni- 
cal X and y axes leaving the Hamiltonian unchanged only 
if the hopping is isotropic. This invariance allows us to 
interchange and ky on the left-hand side of Eq. ([6]) 
resulting in the conditions 

ky = —kx — q + 2m: (7) 

and by reflection symmetry 

— ky = —kx — q + 2mr, (8) 
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where Q = {q,q). For nearest neighbour hopping, the 
resultant condition, k^. ± ky = —it + 2n7r, is the solu- 
tion to cos kx + cos ky = 0, which defines the Fermi sur- 
face for the non-interacting system. If only next-nearest 
neighbour hopping is present, the value of the wave vec- 
tor that leaves the kinetic energy term unchanged after 
a particle-hole transformation is Q = (7r,0)) or (0, tt). 
Coupled with Eq. ^ and reflection symmetry we also 
obtain the Fermi surface of the non-interacting system. 
The interactions need not be the local on-site repulsion 
in the Hubbard model for the surface of zeros to per- 
sist. Nearest- neighbour interactions of strength V de- 
pend only on the particle density and hence are inde- 
pendent of Q under a particle-hole transformation. Such 
interactions renormalise the chemical potential from U /2 
(on-site interactions) to {U + 2V) /2 at the particle- hole 
symmetric point. 

However, an implicit assumption in our proof which 
allows for the interchange of the momenta in Eq. ([5]) 
is that the hopping is isotropic. Nonetheless, the result 
we have obtained is independent of the isotropy of the 
hopping. That is, our proof applies equally when the 
band structure is of the form t{k) = tx cos k^ + ty cos ky, 
with tx ty To prove this, we consider the moments 



M:{k) 



duj 
2^ 



(9) 



of the Green function. For simplicity, we have set h — I. 
Using the Heisenberg equations of motion, we reducei^ 
the moments in real space 

M:{i,j) = I [H--- [H, C.,] • • -Intimcs, cjj) 

+ ({Q„[---[c],,i?]---i/],i/]ntimcJ>] (10) 

to a string of commutators of the electron creation 
or annihilation operators with the Hubbard Hamilto- 
nian. The right-hand side of this expression is evalu- 
ated at equal times. To evaluate the string of commu- 
tators, it sufhces to focus on the properties of — 

[■ ■ ■ [cia,H],- ■ ■ i/]n times, whcrc by construction, K^"^ = 
Cia ■ We write the Hubbard Hamiltonian as H = Ht + Hu 
where Hu includes the interaction as well as the chemical 
potential terms and Ht the hopping term. The form of 
the first commutator. 



suggests that we seek a solution of the form 



K. 



(n) 



Q 



in) 



(11) 



(12) 



where Q-"'' = [• • • [cia,Hu], ■ ■ ■ -ff^ln times involves a string 
containing Hu n times and in Kja , Ht appears at least 
once. Our proof hinges on the form of Q^^^ which we 



write in general as Q^"-* = anCiaiii-^ -\- PnCia. The solu- 
tion for the coefficients 



a„+i = {U - ^J,)a„ + U{-^j,y 



(13) 



-,("+1) 



is determined from the recursion relationship Q\ 

[qI^"* , Hu] . In the moments, the quantity which appears 
is 

{{Qt\c]^}) = % K + /?„] ^ (5.,7„. (14) 
Consequently, the moments simplify to 

M,:(z,j) = %7n + ^5]t.; (({A(:\4})-f /..c.) .(15) 

The criterion for the zeros of the Green function now 
reduces to a condition on the parity of the right-hand side 
ofEq. (jlSp. Consider the case of half-filling, particle- hole 
symmetry and nearest-neighbour hopping. Under these 
conditions, {rii^) — 1/2 and by particle-hole symmetry, 
/i = U /2. The expressions for a„ and /3„ lays plain that 
the resultant coefficients 



In 



C/\" 1 + (-1)" 



(16) 



vanish for n odd. Consequently, Ga{k,uj) is an even 
function if the second term in Eq. (|15p vanishes. In 
Fourier space, the second term is proportional to the non- 
interacting band structure t{k). The momenta at which 
t{k) = define the Fermi surface of the non-interacting 
system. Note, the condition t{k) — which defines the 
surface of zeros is independent of the anisotropy of the 
hopping. We conclude that when particle-hole symmetry 
is present, G{p = pr, 0) = for a Mott insulator, where 
Pi? is the Fermi surface for the non- interacting system. In 
this case, the volume of the surface of zeros is identically 
equal to the particle density. This constitutes one of the 
few exact results for Mott insulators that is independent 
of spatial dimension or at least as long as d ^ oo. As 
mentioned previously, in d = oo, there is no Luttinger 
surface as E diverges as l/to for all momenta^^. Finally, 
the only condition for the applicability of our proof is 
that the form of the spectral function leads to the conti- 
nuity of i?(T(k, Lu) ai uj — 0. Hence, the minimal condition 
is that the spectral function is be continuous at a; = 0. 
Therefore, if there is a gapless quasiparticle excitation, 
for example, Ao-(k, e) = S{uj), our proof becomes invalid. 



B. Away From Particle-Hole Symmetry 

What happens when particle-hole symmetry is broken? 
To consider this regime, we write the electron density 



^2^y hm / — G(fc, 
Vt^o+y 2tt ^ 



loj e 



(17) 
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as an integral of the time-ordered Green function where 
the factor of two counts up and down spin electrons. In 
proving Luttinger's theorem, one uses the identity 

G{k,iLo) = — logG-i(fc,iw) + G'(fc,iu;)— S(/c,iw)(I8) 
aiuj oiuj 

Implicit in this expression is the Dyson equation, 

G"i(fc, iuj) = Go \fc, iuj) + iuj), (19) 

where S is the self energy and Gq the Green function for 
the non-interacting problem. The density, n — Ii — I2, is 
now a sum of two terms 

/OC 7 o 
^^-logG-(M.) (20) 

/OO 7 o 
^G(fc,^c.)^E(fc,^c.). (21) 
-00 27r diLu 

Luttinger— proved that I2 vanished and hence the elec- 
tron density is given simply by Eq. H]). In fact, 
Dzyaloshinskii^ has claimed that I2 vanishes for a Mott 
insulator. Central to this proof is the existence of a per- 
turbative expansion for the self energy around the atomic 
limit. Based on the self-energy, the Luttinger- Ward (LW) 
functional^, 



(22) 



can be constructed which for a Fermi liquid has a pertur- 
bative expansion in terms of skeleton diagrams. In gen- 
eral, the LW functional is assumed to have a perturbative 
expansion. As such any perturbative LW functional must 
be free of singularities and vanish as uj 00. We show 
here that for the Mott problem, singularities arise and 
it is precisely from the singularities that a breakdown of 
Luttinger's theorem arises. 

To see how Luttinger's theorem fails for a Mott insu- 
lator, consider the exact temperature-dependent Green 
function, 



G{k,iuj) 



iuj + fi + U/2 — 'Eioc{iuj) 

iiu + fi 
(iuj — Ei)(iu! — E2)' 



in the atomic limit where 



Sloe = 17 + 17 , 

2 \ 2 J IUJ 



1 



(23) 



(24) 



and Ei^2 = — ^± C//2. As the Mott gap is well-formed in 
this limit, any conclusion we reach regarding I2 will hold 
as long as U ^ t. Can Eq. be constructed from 

the non-interacting limit? To all orders in perturbation 
theory at T = 0, the self-energy is given by 



Such a self-energy cannot describe the two-peak struc- 
ture of the Mott insulating state. That is, starting from 
the non-interacting system, one cannot obtain the Mott 
gap perturbatively. It is this breakdown of perturbation 
theory in generating the Mott gap that is central to the 
ultimate breakdown of the Luttinger sum-rule on the vol- 
ume of the surface of zeros. Given that Spert 7^ Sioc, the 
corresponding LW functional, ^{G) cannot be obtained 
perturbatively. Consequently, we must resort to a non- 
perturbative method to construct the LW functional. To 
gain some insight into what the corresponding LW func- 
tional looks like, we rewrite Sioc in terms of the exact 
G, 



l±Vl + U^G^ 
2G 



(26) 



by eliminating the 1 / {iLUn + m) factor by using the Dyson 
equation and solving the subsequent quadratic equation. 
Here, the upper and lower signs should be used when 
> U/2 and luj+fi] < U/2 respectively. As i;[G] sat- 
isfies the Dyson equation, it is exactly given by the func- 
tional derivative of the exact LW functional with respect 
to G. Hence, Eq. (gSD implies that we know 5<^[G]/5G 
at the saddle point of $. Constructing $ in general, 
however, requires complete knowledge of 5^[G]/6G not 
simply at the saddle point. How then do wc construct 

For the problem at hand, there are two requirements 
that any approximate expression for ^[G] must satisfy: 
1) it must contain a singular part and 2) I2 computed 
from any approximate LW functional must agree with 
a direct calculation based on the second of Eqs. (|20p . 
The approximate LW-functional we derive here satisfies 
both of these requirements and hence lends credence to 
the method. Faced with the similar problem of obtaining 
^[G] knowing only its exact derivative at one particular 
value of G for a spin-density wave problem, Altshuler, 
et aliii simply integrated S[G] with respect to G to ob- 
tain an approximate LW functional. They showed that 
this procedure to be internally consistent for their prob- 
lem as I2 evaluated with the approximate LW functional 
agreed with a direct calculation of I2 from Eq. (PO)) . We 
adopt this approach and check its internal consistency in 
a similar manner. The integral of 5][G], 



$(iw) 



^1 ^l + U^G^{iLo) 
±7: log — , 

2 \.JTTTPGHiu) + 1 

$rcg(it^) + $smg(iw). 



I 



contains both a regular as well as a singular part. 



$sing(it^) 



1 G^jiuj) 

2 °^ G[iuo) ' 



(27) 



(28) 



-'pert 



u. 



(25) 



Although the atomic limit differs from the spin-density 
wave problem treated by Altshuler, et alJ^, the approx- 
imate LW functionals are identical. This state of affairs 
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obtains because of the similarity between the self ener- 
gies of the two problems. To evaluate I2, we note that 
only the singular part of ^[G] contributes. The result 



-2i 



27r 



2e(/i)-e( 



diuj 

-Ei)- 



(29) 



is in agreement with a direct calculation^*^ of I2 based on 
the second equation in Eq. (^0)1 . This agreement suggests 
that our approximate expression for <&[G] captures the 
essence of the breakdown of Luttinger's theorem. Note 
I2 term vanishes only in the presence of particle-hole sym- 
metry (/i = 0) . The modified Luttinger theorem becomes 



(Pk 

'(o,fc)>o (27i-)2 



(30) 



The modified Luttinger theorem is also valid even in the 
presence of small hopping, that is, U t, as can be seen 
by considering 



1 



iLO — t{k) + /i + C//2 — Tiioc{iijj) 

iuj + ^ 
{iuj - Ei{k)){iuj - E2{k)) 



(31) 



for the Green function. Here t{k) is the Fourier trans- 
formation of the hopping element, tij. Substitution of G 
and Yj into the second term of Eq. PT|) leads to 



h = 



2e(^) - j 



<pk 



[Q{-E^{k)) + Q{-E2{k))]. (32) 



Once again, I2 vanishes in the presence of particle-hole 
symmetry and hence can be rewritten as Eq. (j29p . 

The form of the singular part of the Luttinger- Ward 
functional compels a simpler formulation of the electron 
density. Since the regular part of I2 vanishes, we can use 
Eqs. (PO)) and ([^5]) to recast the electron density as a 
sum of two contributions, 

n ^ J ^-^logG{k,iLu)Go{k,iuj) 



(fk 



(fk 



G{k,uj=0)>0 (271")^ J Gt=o,u=o{k,uj=0)>0 {'^''^Y 

(33) 

each of the Luttinger form. Although both of the terms 
in this expression (as well as in Eq. (30)) contain discon- 
tiuities, the discontinuities cancel in the sum leading to 
the density being a continuous function of the chemical 
potential. Eq. (|33p is valid in the atomic limit as well as 
in the small hopping regime and represents the general 
form of Luttinger's theorem for a Mott insulator. Note 
in the weak hopping limit, only the first term differs from 
that in the atomic limit. In interpreting Eq. (|33p . it is 
important to remember that the second term is not equiv- 
alent to I2 ■ Part of I2 cancels one of the Green functions 



in the first term. The term which is left over accounts for 
the fact that the chemical potential can be placed arbi- 
trarily within the gap as emphasized previouslyi^. This 
ambiguity, of course, is absent for a soft gap as in the case 
of the pseudogap in the doped case. In this case, how- 
ever, the exact self-energy is not known and no recasting 
of the Luttinger theorem as the general statement in Eq. 
([55)) is possible. What the current analysis lays plain is 
that the singular part of the LW functional, which is ab- 
sent for a Fermi liquid, leads to the break-down of the 
Luttinger sum rule on the surface of zeros in the absence 
of particle-hole symmetry. 

The current analysis can be extended to finite temper- 
ature. At finite T, I2, 



h{T > 0) 



+/(-M 



-2/(m),(34) 



can be evaluated using Matsubara frequencies and the 
singularity which originally existed on the real frequency 
axis can be removed. Eq. goes over smoothly to the 
zero-tcmperature l2{T = 0) evaluated by real frequency 
integration. Therefore, the non- vanishing of I2{T — 0) is 
independent of the regularization of the singular part of 
the LW functional and is a generic feature of a Mott in- 
sulator. This result is significant because Dzyaloshinskii^ 
proposed that the non- vanishing of I2 stemmed from the 
method Altshuler, et. alf^i used to regularise the singu- 
lar integrals. Namely that the breakdown of Luttinger's 
theorem for a spin-density wave arises from a T = regu- 
larisation of divergent integrals that can only be reached 
if there is a T = phase transition. The current work 
establishes that there is no such phase transition and the 
T = result is connected adiabatically to the finitc-T 
result. Gonsequently, the breakdown of Luttinger's theo- 
rem lies elsewhere. In both the Mott insulator and spin- 
density wave problems, no perturbative expansion ex- 
ists for the self-energy around the non- interacting limit. 
Should this fail, there will always be a singular part of the 
LW functional and I2 will be finite. Altshuler et. alJ^ 
made a connection between such a break down and the 
chiral anomaly in particle physics. While at the atomic 
limit, this association might be appropriate, it is unclear 
whether this analogy holds for the general case. To reit- 
erate, for a Mott system lacking particle-hole symmetry 
but possessing a divergent self-energy, the singular part 
of GdujY will always integrate to a non-zero value. 



C. Consequences 

Several claims follow necessarily from these results. 
CI. There are no non-trivial zeros of the single-particle 
Green function in the single-impurity Anderson model. 
Because the gap is replaced by the Kondo resonance, no 
zeros of the Green function obtain for the single- impurity 
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problem. 

C2. At the surface of zeros, the self-energy at zero fre- 
quency diverges. Write the single-particle time-ordered 
Green function as Ga{k, oj) — 1/{lu — e{k) — ^T,„{k, uj) — 
i^'Sa-{k,uj)), where S is the self-energy. For a particle- 
hole symmetric band structure, the single-particle Green 
function vanishes linearly at the Luttinger surface, k^: 
Ga{k,uj = 0) = const, x (fc — fci). This implies that 

5I](fc,0) oc (5(fc- fci). (35) 

Note, however, that 5G(fc,u;) = for all energies within 
the gap. By inverting the Green function, it follows that 

5RE<,(fc,w) cx (fc- (36) 

proving C2. To reiterate, in d = oc^^ no Luttinger sur- 
face exists as the self-energy diverges for all momenta at 
Lu = 0. While such a divergence is appropriate for d = oo, 
it is clearly unphysical for a finite-dimensional system. 
The divergence in Eq. (|36p prevents the renormalized 
energy band E{k) — e(fc) -I- ^EaikjUi) from crossing the 
Fermi energy. The result is an insulating state. Indeed, 
in numerical studies^^ on the Hubbard model at half- 
filling with nearest-neighbour hopping, Scr(7r,0,[i; — 0) 
has been observed to diverge as our theorem indicates 
it must. However, Jarrell and co-worker^i^ attributed 
antiferromagnetism as the cause of the divergence. Our 
theorem indicates that the zeros are independent of the 
ground state (be it ordered or not as in the case of a spin 
liquid) as long as the Mott gap is present. The zeros are 
a direct consequence of Mottness itself. That Mottness 
and zeros are one and the same indicates that the di- 
vergence of the self energy provides a general mechanism 
for insulating states in the absence of broken symmetry. 
That is, in a finite-dimensional lattice, the divergence 
of the self-energy at a contiuously connected momentum 
surface is the general mechanism by which the Mott in- 
sulating state obtains. 

C3. Zeros represent a breakdown of weak-coupling per- 
turbation theory. This follows directly from C2. A di- 
vergence of the self-energy is the general signature of the 
breakdown of perturbation theory. Zeros offer a concrete 
way of realising this breakdown. In d = 1, this break- 
down occurs for all U ^0. In d = 2 in the particle-hole 
symmetric case, the critical value of U is not known, 
though all numericsi^ii^ indicates that the only special 
point is U = 0. 

C4. The surface of zeros of the single-particle Green 
function is absent from projected models at half-filling. 
Since it is common in the study of doped Mott insula- 
tors to use projected models, it instructive to evaluate 
whether such truncations admit a surface of zeros. Pro- 
jecting out double occupancy, as in the t—J model, erases 
the spectral weight above the chemical potential at half- 
filling. Consequently, the real part of the Green function 
reduces to the first integral in Eq. which is always 
non-zero. C4 is thus proven. Transforming the operators 
in the t—J model to respect the no double occupancy 



condition is of no help as the problem stems from the loss 
of spectral weight above the gap once projection occurs. 
As the surface of zeros occurs at zero energy and is lo- 
cated in momentum space, it should certainly be present 
in a low-energy theory of the Hubbard model. However, 
it is clear from Eq. ^ that zeros of the Green function 
stem from a sum rule^^ connecting low and high ener- 
gies. Hence, it is a priori expected that the zero surface 
would be sensitive to the retention of the spectral weight 
at high energies. 

The absence of zeros in the t—J model at half-filling 
is in actuality related to the problem of the robustness of 
zeros and the location of the chemical potential at half- 
filling in a Mott insulator. At T = 0, the chemical po- 
tential is a free parameter that can be located anywhere 
in the Mott gap. Consider, the extreme case of placing 
the chemical potential atop the lower Hubbard band and 
sending U to infinity. In this case, the integrand in the 
second term in Eq. (4) has an infinite energy denomina- 
tor and hence the second term vanishes. Consequently, 
there are no zeros in this case. The actual realisation of 
this is the t—J model at "half-filling."— Hence, there are 
certain locations for the chemical potential for which the 
zero line vanishes. This does not diminish the significance 
of the zero line as the defining feature of a Mott insula- 
tor, however. That the chemical potential is arbitrary at 
T — Q indicates that the T = value of the chemical po- 
tential is not a defining feature of a Mott insulator. What 
is the defining feature of a Mott insulator is that at half- 
filling, half the spectral weight lies above the gap. Such 
a schism in the spectral weight guarantees that the real 
part of the Green function must change sign along some 
momentum surface for some energy or range of energies 
within the gap as emphasized by DzyaloshinskiS. Our 
claim that the surface of zeros defines the Mott insulator 
is simply that dynamical generation of a gap, which at 
half filling results in half the spectral weight lying above 
and below the gap, leads to a sign change of SRG for some 
(not necessarily all) energies within the gap. In this vein, 
the t—3 model is not a realistic model of the Mott state 
because a zero line is strictly absent. 
C5. Even at infinitesimal doping, the t-J and Hubbard 
models probably do not yield equivalent values for the 
quasiparticle weight. Because the chemical potential sits 
atop the lower Hubbard band in the t—J model at "half- 
filling"^^, perhaps the proper way to compare the with 
the Hubbard model is in the limit of infinitesimal dop- 
ing (see Figs. (|3^) and (c)). Numerical and analytical 
studies on the one hole syste m^°i^^i^^ find a quasiparticle 
in the t-J model with weight J/t at (7r/2,7r/2) whereas 
in the Hubbard model^^, the quasiparticle weight van- 
ishes as Z (X , 9 > 0, L the system size. Variational 
calculations^i also yield a finite value of Z in the extrap- 
olated limit of n = 1^. While none of this constitutes 
a proof, it is highly suggestive that the value of Z is 
tied to the presence of the upper Hubbard band as has 
been emphasized previousl y^^i^^ . In the one-hole system, 
sufficient spectral weight must lie above the chemical po- 
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tential for Eq. (4) to vanish. There is no guarantee that 
this state of affairs obtains for the t — J model since no 
spectral weight was above the gap at the outset. No such 
problem arises for the Hubbard model. 
C6. // a Mott gap opens, zeros of the single-particle 
Green function still persist when the particle-hole sym- 
metry is broken weakly. As remarked earlier, all that 
is necessary to establish is that for some energies within 
the gap, the real part of the Green function changes sign. 
At present, our proof applies to any kind of band struc- 
ture that is generated from hopping processes which re- 
main unchanged after the application of Eq. In 
general, the two kinds of hopping processes transform as 
e(7r - kx.TT - ky) = -e{kx,ky) and i'(7r - kx,n - ky) = 
t'{kx,ky). The latter describes next-nearest neighbour 
hopping and as is present in the cuprates. If only such 
hopping is present, the surface of zeros is no longer the di- 
agonal (tt, 0) to (0, tt) (or the point 7r/2 in ID), but rather 
the "cross" (0,7r/2) to (7r,7r/2) and in/2,0) to (7r/2,7r) 
(or, in ID, the points —tt/A and 37r/4). When both types 
of hopping are present, no symmetry arguments can be 
made. Our proof in this case will rely on a key assump- 
tion: the Green function is a continuous function of the 
hopping parameters t and t' . Hence, strictly speaking our 
proof applies only when i' <C t. When only t is present, 
Ra{k,0) has one sign (plus) near k = (0,0) (or, in ID, 
k=G), and the opposite (minus) near k ~ (7r,7r) (k — ix 
in ID) and will vanish on the zero line. Alternatively, 
if we have t' hopping, R^{k,0) will have a certain sign 
near k = (0,0) and k — (7r,7r), and the opposite sign 
near k = (0, tt) and k — (tt, 0) and will vanish on the 
"cross". From continuity, for t' <C t, Rc{k;t,t') will have 
the same sign structure as Rcr{k;t,t' — 0). That is, it 
will change sign when going from (0, 0) to (tt, tt) regard- 
less of the path taken. Therefore, the line of zeros exists 
for small enough t', the relevant limit for the cuprates. 
In the opposite limit, i' 3> a similar argument holds. 
Whether a proof exists for the case of a strong violation 
of particle-hole symmetry is not known. 



III. UTILITY OF ZEROS: PSEUDOGAP PHASE 
OF DOPED MOTT INSULATORS 

Ultimately, the utility of the surface of zeros will be de- 
termined by its experimental relevance. As we have seen 
above, the volume of the surface of zeros can only be cal- 
culated explicitly in the case of particle-hole symmetry. 
At finite doping where this symmetry is explicitly bro- 
ken, the existence (though not the volume) of the zero 
surface can nonetheless be established. Two indepen- 
dent arguments are relevant here. The first is based on 
the distribution of spectral weight in the spectral func- 
tion and the other on the fact that Fermi arcs, as seen 
experimentall y^^i^^ in the doped cuprates, necessitate the 
existence of a surface of zeros. The general arguments 
made here for the interdependence of zeros and Fermi 
arcs augment the numerical evidence found by Stanescu 



and Kotliar— for the same effect. 

To illustrate that the spectral weight distribution in 
a lightly doped Mott insulator supports a zero surface, 
we consider the spectral function shown in Fig. ([T]). The 
computational scheme used to produce this spectral func- 
tion has been detailed elsewhereiS. and is in agreement 
with results from state-of-the-artJ^ calculations on the 
2D Hubbard model. Two features are relevant. First, 
it possesses a depressed density of states at the chemi- 
cal potential for a wide range of momenta. This leads 
to a density of states which vanishes algebraically at the 
chemical potential, as is seen experimentally^^. Such a 
dynamically generated pseudogap which occurs without 
any symmetry breaking has been confirmed by all re- 
cent numerical computations on the the doped Hubbard 
model^iiii^i^i^a. Hence, that 9G(0,p) = along some 
contour in momentum space for a doped Mott insula- 
tor is not in dispute neither theoretically nor experimen- 
tally What about 3fiG'(0, p)? As is clear from Eq. g]), 
5RG(0, p) = if along some contour in momentum space, 
the spectral weight changes from being predominantly 
below the chemical potential to lying above. At half- 
filling, the zero surface obtains entirely because most of 
the spectral weight at (tt, tt) lies above the chemical po- 
tential, whereas at (0, 0), it lies below. As is evident, this 
trend still persists for x « 1 as Fig. ^ attests. Hence, 
Ra{k,0) still has the same sign structure as in the un- 
doped case. Consequently, a zero surface must exist. 

Ultimately, satisfying the zero condition, Eq. re- 
quires spectral weight to lie immediately above the chem- 
ical potential. Spectral weight transfer'^^ across the Mott 
gap is the mediator. The weight of the peak above the 
chemical potential scales as 2x -\- f{x, t/U)^ (strictly 2x 
in the t-J model) while the weight below the chemical 
potential is determined by the filling, 1 — a;. Whether 
or not the redistributed spectral weight is symmetric or 
not around the chemical potential will determine how 
severely the Lutinger volume is violated. There are only 
two options as depicted in Fig. 

Weak violation of Luttinger volume: In order to 
satisfy Luttinger theorem, the surface of zeros has to be 
close to the (0, tt) — (tt, 0) line (assuming that is small). 
Consider a point on the zero line in the vicinity of (0, tt). 
As is evident from Fig. ([T|) the spectral weight immedi- 
ately above /i (in the vicinity of the (7r,0) — (0,7r) line 
is small compared with the spectral weight below fi and, 
in order for diG to vanish, the chemical potential has to 
be positioned asymmetrically inside the pseudogap (see 
panel A of Fig. 

Strong violation of Luttinger volume: The other 
possibility is that the Luttinger theorem is strongly vio- 
lated and the surface of zeros is somewhere in the vicinity 
of (7r,7r). In that region, the spectral weight of the lower 
Hubbard band is greatly reduced as shown in Fig. 
Consequently, most of the spectral weight lies in the up- 
per Hubbard band and a more symmetrical distribution 
around the chemical potential is possible (see panel B in 
Fig. 0. 
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FIG. 1: Spectral function for a doped Mott insulator at a fill- 
ing of n = 0.95 at T = 0.07t for a path in momentum space 
from (0,0) to (tt.O) to (tt.tt) and then back to (0,0). Eq. dH) 
must undergo a sign change by passing through zero because 
at (0, 0) most of the spectral weight lies below the chemical 
potential, whereas at (7r,7r), it lies above. The temperature 
dependence of the density of states at the Fermi level is shown 
in the inset. That the spectral features are broad near the 
chemical potential is a direct consequence of the divergence 
of the self energy at a; = at the Luttinger surface. The spec- 
tral function was computed using the two-site self-consistent 
method of Stanescu and Phillips^- . "Color Online" 
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To decide between options A and B in Fig. ([2]) we 
appeal to experiments. Two observations support option 
A. First, the pseudogap is in fact asymmetrical^^. Sec- 
ond, consider the recent photoemisison experiments^^ in 
which the temperature dependence of the Fermi arcs has 
been measured. Experimentally, lightly doped cuprates 
possess Fermi arcs^^ along the zone diagonal in the vicin- 
ity of (7r/2,7r/2). Whether the Fermi arcs represent a 
finite T precursor of a Fermi surface and hence quasipar- 
ticles as claimed by some^ can be settled by temperature- 
dependent ARPES experiments in the pseudogap regime. 
KanigelSi et al. performed such temperature-dependent 
ARPES measurements on Bi2Sr2CaCu208+5 (Bi2212) 
and concluded that the Fermi arc length shrinks to zero 
as T/T*, where T* is the temperature at which the 
pseudogap feature appears. Hence, the only remnant 
of the arc at T = 0, is a quasiparticle in the vicinity 
of (7r/2,7r/2). Consequently, Kanigel, et al.— argue for 
a nodal metal. That a nodal metal or a quasiparticle 
band existing over a finite connected region in momen- 
tum space not extending to the zone boundary (as would 
be the case in a T = Fermi arc) cannot be understood 
without the existence of a surface of zeros can be seen 
as follows. Assume a quasiparticle exists at (7r/2,7r/2). 
Then 5RG(0,p) must change sign for all momenta less 
than or greater than {tt/2, 7r/2) as depicted in Fig. ^jP). 
Consider traversing a path through (7r/2,7r/2) and then 
returning along a path that does not cross this point. 
To end up with the correct sign for 5RG, the return path 
must intersect a line across which ^G{0, p) changes sign. 
Since there are no infinities, except at (7r/2,7r/2), the 
only option is for a zero line to exist. The zero line must 
emanate from the (7r/2,7r/2) point and touch the edges 
of the Brillouin zone close to (tt, 0) and (0, tt). A zero 
surface terminating close to (tt, tt) is not an option as 
this would permit the existence of paths that traverse 
the zone diagonal without changing the sign of 3fiG(0, p). 
This would suggest that the zero surface in the doped 
cuprates preserves the Lutinger volume and option A in 
Fig. ([2]) is more consistent with experiment. As a con- 
sequence, Fermi arcs are direct evidence that zeros of 
the single particle Green function must be present in the 
doped cuprates. In a recent paper, Stanescu and Kotliar— 
have argued based on numerics for such an interdepen- 
dence. 



FIG. 2: Schematic depction of two possibilities for the dis- 
tribution of spectral weight in the single-particle density of 
states, n(aj), in a lightly doped Mott insulator. In both cases, 
the peak above the chemical potential represents the low- 
energy spectral weight. The weight of this peak increases 
as least as fast as 2x, where x is the number of holes. A) 
The Luttinger surface lies along the zone diagonal and the 
transferred spectral weight must be asymetrically located rel- 
ative to the chemical potential to lead to a vanishing of 
B) The spectral weight redistribution is symmetrical and the 
Luttinger surface lies away from the zone diagonal, for exam- 
ple close to (tt, tt). 



IV. CONCLUDING REMARKS 

As we have seen the experimental utility of zeros of 
the single-particle Green function is in their relevance 
to ARPES. One of the hallmarks of the normal state of 
the cuprates is an absence^^ of electron-like quasiparti- 
cles. Quasiparticles require a vanishing of the renormal- 
ized band, E{k) = e{k) -\- 3?I](fc,w = 0); but because 
diverges along the surface of zeros, no quasiparticles 
form and broad spectral features are inevitable as seen 
in ARPES in the cuprates^^. The clearest experimental 



9 




n<1 
(metal) 



(0,0) 



(;r,0)(0,0) 



(;r,0) 



FIG. 3: (Color online) Evolution of the surface of zeros in 
the first quadrant of the FBZ. Yellow indicates KG > while 
blue KG < 0. The Hubbard model is constrained to have a 
surface of zeros as n — > 1 whereas the t-J model is not. The 
two options upon doping represent weak-coupling (la or c to 
d) and strong-coupling (la to b). The transition from (a) 
to (b) requires a critical point at ric whereas (d) does not. 
Experiments^!^ indicating an insulating state for n > ric are 
consistent with an abrupt transition from (a) to (b). 



signature that the surface of zeros exists is the recent 
temperature- dependent ARPES experiments that indi- 
cate that the Fermi arcs shrink to a point as T — > 0. Since 
the surface clearly exists for the cuprates, the only out- 



standing question is how does the surface of zeros evolve 
as a function of doping. Various options are shown in 
Fig. An abrupt transition from a surface of zeros 

to a Fermi surface would describe a transition from an 
insulator to a metal. Such a transition would require 
a phase transition at Xc, the doping level at which the 
pseudogap terminates. Alternatively, quasiparticles and 
zeros could co-exist. While we have advocated the for- 
mer scenario based on a calculation of the conductivity 
which reveals that the pseudogap is an insulating gap'^^, 
a result consistent with experimen1*^i2i ultimately, both 
scenarios are possible, in principle^. The former cor- 
responds to an insulator (or nodal metal) whereas the 
latter describes a metal. The recent angle-resolved pho- 
toemisison experiments^ seem to indicate that the only 
co-existence of quasiparticles and zeros occurs at a single 
point indicating that a Fermi surface is possible only for 
some doping level exceeding depicted in the upper 

panels in Fig. ([3]). 
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